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Abstract 

The Z^-critical parabolic-elliptic Patlak-Keller-Segel system is a classical model of chemo- 
tactic aggregation in micro-organisms well-known to have critical mass phenomena [S] . In 
this paper we study this critical mass phenomenon in the context of Patlak-Keller-Segel models 
with spatially varying diffusivity of the chemo-attractant. The primary issue is how, if possi- 
ble, one localizes the presence of the inhomogeneity in the nonlocal term. We also provide new 
blow-up results for critical homogeneous problems with nonlinear diffusion, showing that there 
exist blow-up solutions with arbitrarily large (positive) initial free energy. For several non-trivial 
technical reasons, which we discuss in more detail below, we work in dimensions three and higher 
where L 1 -critical variants of the PKS have porous media-type nonlinear diffusion on the organism 
density, resulting in finite speed of propagation and simplified functional inequalities. 



1 Introduction 

The most widely studied mathematical models of nonlocal aggregation phenomena are the parabolic- 
elliptic Patlak-Keller-Segel (PKS) models, originally introduced to study the chemotaxis of microor- 
ganisms |37[ [28l [26l [25] ■ In this paper we consider L 1 -critical cases of the form, 

u t + V • (uVc) = Au 2 ~ 2/d 

-V • (a(x)Vc) + 7(x)c = u (1.1) 
u(0,x) = u (x) G L 1 + (R d ;(l + \x\ 2 )dx)nL°°(R d ), d>3, 

where L+(M. d ; n) := {/ G L 1 (R ; //):/> 0}. We define m := 2— 2/d which is the nonlinear diffusion 
power for which the system is L 1 -critical. In this context, L 1 -critical refers to the approximate 
balance of the opposing forces of diffusion and aggregation in the limit of L 1 concentration, indicating 
that there must be a non-zero, but finite, amount of mass concentration at any possible blow-up. 



As weak solutions to (1.1) conserve mass, we will henceforth refer to ||uo||i = ||u(£)||i = M{u). For 



all our work, we assume that a(x) is strictly positive which ensures the PDE for c to be uniformly 



elliptic. The model (1.1) is a generalization of the classical parabolic-elliptic 2D PKS model which 
has received considerable attention over the years (see the review [26J and [271 [22], [Til QUI IE])- 
For aggregation occurring in three dimensional environments, the nonlinear diffusion models an 
over-crowding phenomenon which retains the L 1 criticality, unlike the corresponding supercritical 
model with linear diffusion. It is well-known that all such L 1 critical models exhibit critical mass 
phenomena: there exists some M c > such that if M(u) < M c then every solution exists globally, 
and if M(u) > M c then there exist solutions which blow up in finite time (see for instance |36} ITTl 
[TOl l9l [2]). We refer to the special case of a{x) = 1, j(x) = as the scale-invariant problem, because 
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solutions are invariant under an L 1 scaling in space and this scaling symmetry plays a fundamental 
role in the global theory. 

In this paper we estimate the critical mass, and under certain restrictions, show that this esti- 



mate is sharp for the PKS model (1.1) with spatially variable coefficients in the chemo-attractant 
PDE. The study of spatially variable coefficients raises questions which are both mathematically 
interesting and relevant for biological applications where chemotaxis occurs in a spatially inhomo- 
geneous medium. We introduce inhomogeneity specifically to the nonlocal nonlinearity, since it is 
not obvious that the effect on the qualitative behavior should be localized. One could also con- 
sider inhomogeneous mobility for the transport of u, but this does not add any new interesting 
complications, so for the sake of clarity we leave this generalization out (see Remark [7| 

For many problems which are scale-invariant, determining critical thresholds is straightforward 
with a classical procedure: sharp functional inequalities to prove global existence below the threshold 
and virial methods to prove blow-up above (see more discussion below). As mentioned above, 
inhomogeneous critical problems can introduce new challenges and change the qualitative behavior 
of models. The inhomogeneity will alter the behavior of the solution at all length-scales differently 
and well-established methods for homogeneous problems can potentially break down. When proving 
global existence, the idea is to rely on an approximate scale-invariance which is recovered in the limit 
of infinitesimal length-scales. Naturally, the challenge in this kind of localization lies in the control 
of the nonlocal error. On the other hand, this viewpoint does not help at all when studying blow-up 
dynamics, as one needs to exploit some monotonicity which persists as the solution blows up and is 
valid for nonlocal, macroscopic interactions. Most approximations, including the available results 
involving mass comparison, will introduce unbounded errors as the solution blows up, making it 
challenging or impossible to apply straightforward monotonicity estimates. This issue makes the 
blow-up problem above critical mass significantly more subtle than the problem of global existence 
below critical mass. 

For the inhomogeneous models, our results state that blow-up is spatially localized: it only 
depends on the local (minimum) value of a(x). We estimate the critical mass to be given by 



M c 



(m - l)C*Cd 



where C* is the optimal constant in the Hardy-Littlewood-Sobolev inequality discussed in |9j (see 
( |1.6[) be low) and Cd is the normalization constant in the Newtonian potential, given explicitly below 
m ~JTI\). We prove that if M(u) < M c then the solution is global and uniformly bounded in 
L°°((0,oo) x M. d ) (Theorem Further, we prove a stronger blow-up criteria which is also valid 
in supercritical mass regimes and shows that the critical mass is truly dependent on space, which 
accurately confirms the localized nature of the blow-up (Theorem [6]). More precisely, we show that 
in order for a solution to blow-up at a time T* £ (0, oo] , for every sequence of times tk there exists 
a subsequence (not relabeled) and a sequence of points Xk such that 



f ( 2a(x k ) \ d/2 
limsuplimsup / u(tj c ,x)dx > liminf I - — — — I 

r _> + A-^oo J\x-x k \<r k ^°° V( m- l)k* c <i/ 



Combining this with mass comparison methods, we prove that radially symmetric solutions with 
M(u) = M c are global and uniformly bounded. Note, the uniform bound is in contrast to the 
infinite time aggregation exhibited in IR 2 [TO] (see more discussion below) . In the direction of finite 
time blow-up we prove that if 7 = and a(x) is radially symmetric and monotone non-decreasing 
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in a neighborhood of the origin, then for all M > M c we may construct a solution with M(u) = M 
which blows up in finite time (Theorem [§]). See the discussion below for the proof of the blow-up 
result and the relevant construction of barriers. 



A key quantity in the study of (1.1) is the dissipated free energy, given by 



F(u) = [ u m (x)dx - 1 / u(x)c(x)dx. (1.2) 

m- 1 J 2 J 

The first term is usually referred to as the entropy and the latter term is referred to as the inter- 



action energy or potential energy. Formally, (1.1) is a gradient flow with respect to the Euclidean 



Wasserstein distance for ( |l.2[ ) (see e.g. [7]), but this will not be relevant for our work (however, see 
|8j for work on the threshold problem with linear diffusion where this structure is the key tool). 

We consider the question of global existence first. It is well-known that solutions to systems 



such as (1.1) exist as long as they remain equi-integrable (27J HE! EJ [2] . The use of sharp functional 
inequalities to identify when a mixed-sign energy such as ( |1.2[ ) is coercive is classical, for example 
|43| ITT| [9j |2]. In the context of classical PKS and similar models, sharp Hardy-Littlewood- 
Sobolev inequalities prove that when M{u) < M c the free energy uniformly controls the entropy 
(in this case the L m norm), which rules out any loss of equi-integrability (see [221 ED El El 12] ) • In 
[2], it was shown for relevant systems with c = K, * u, for interaction kernels K, in a general class, 
that the critical mass is governed only by the asymptotic expansion of /C at the origin. There, 
the asymptotic expansion describes the singular contribution to the potential energy J ucdx and 
the remaining error is 'subcritical' in the sense that it can be controlled by norms weaker than 
the entropy. This is, in spirit, the approach taken here to prove Theorems [5] and [6] in |2] however 



a different method must be applied for (1.1), since in our case c is not given by a convolution. 
Instead, we use a pseudo-differential operator to explicitly evaluate the leading order contribution 
to the potential energy as a quadratic form and show that the remainder is subcritical. This allows 
to refine the classical methods and estimate the critical mass. Additionally, the critical mass blow- 
up criterion (Theorem [6]) is deduced by refining and iterating the arguments of |10| (see below for 
more information). A crucial new tool here is a geometric covering lemma which, combined with 
the concentration compactness principle [32J, allows to decompose blow-up solutions into isolated 
regions in which a localized energy argument can be applied. 



Conventionally, the simplest way of proving blow-up for systems such as (1.1) above the critical 
threshold is the well-known virial method, used in, for example |36} HOI El EE]- Applied to the 
scale-invariant problem in d > 3, this method consists of two steps (see e.g. j9]). First, one 
uses supercritical mass and a scaling argument to construct initial data with negative free energy. 
Second, one shows that negative free energy would force the second moment to zero in finite time. 



Although we do not claim that this cannot be done for systems such as (1.1), approximating the 
nonlocal term with pseudo-differential operators as done in the global existence results introduces 
error terms for which there seems to be no obvious way to control into blow-up. Hence, we instead 
use a different finite time blow-up proof which is able to treat the local and nonlocal properties of 
the solution with sufficient precision. 

The proof of finite time blow-up is based on a maximum-principle type argument on the mass 



distribution of solutions. The use of mass comparison principles for models such as (1.1) is by now 
classical [HI [201 [381 123 EQl [16]. However, our work is the first one (to our knowledge) which 
treats blow-up phenomena with mass arbitrarily close to the critical threshold. This is possible 



due to the precision of the barrier. Specifically, we compare solutions of (1.1) to the extremizers 
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of the sharp Hardy-Littlewood-Sobolev inequality which governs the critical mass, which are also 
stationary solutions of the scale-invariant PDE (Theorem [4]) . More details are discussed in Section 

Due to the nature of the near-critical regime we address, the mass comparison argument we adopt 
is rather delicate and depends on certain regularity properties of solutions as well as the barrier. 



Among other complications, the degenerate diffusion in (1.1) implies that classical regularity is not 
available everywhere and the support of solutions move at a finite speed. In order to deal with 
the generic presence of the free boundary, we utilize the viscosity solution theory developed for 
degenerate diffusion equations with drift (see [29]) to prove that strictly positive solutions remain 
positive until blow-up (see Appendix). Then we need to make a careful approximation argument 
with smooth, strictly positive solutions, as in |3U| . We add that here extra care must be taken due 
to the finite time blow-up (Lemma [7]). 

With some modification, this blow-up proof can also provide estimates from above on how quickly 
mass concentrates, providing also a lower bound on the blow-up time by use of a mass supersolution, 
which is perhaps a more common use of mass comparison principles. Indeed, this is how the global 
existence at critical mass is proved for radial solutions (see below) . More interestingly, when applied 
to the scale-invariant problem, our method yields blow-up for a class of initial data that the results 
in [9], obtained via a virial method, do not cover (see Theorem [9] below) . In particular, we exhibit 
radially symmetric solutions with initially positive free energy that concentrate in finite time. This 
seems to indicate that the approach we employ could potentially provide different kinds of blow-up 
results elsewhere, even when straightforward virial methods can be applied. 

In the last section of the paper we prove that under the assumption of radial symmetry, solutions 
with critical mass M(u) = M c exist globally and are uniformly bounded. Unlike the homogeneous 
case, one cannot expect radial monotonicity to hold (even when a(x) is radial), which means that 
the construction of a mass supersolution is not enough to deduce such a result without additional 
information. On the other hand, since the solutions have exactly critical mass, a simple corollary 
of Theorem [6] below is that in order to blow up, such threshold solutions must concentrate all of 
the mass into a single point where the minimum of a(x) is achieved. Note that Theorem [6] does not 
imply such rigid behavior at blow-up for solutions with larger than critical mass. Mass comparison 
methods similar to those employed to prove finite time blow-up and the arguments in j5j 1431 1381 [30] 
are used to show that this concentration cannot occur in finite or infinite time. Hence, the solution 
must be global and uniformly bounded. 

We restrict ourselves to the case M. d for d > 3 for several reasons. Firstly, note that the d > 3 
case is qualitatively different than the d = 2 cases aside from the finite speed of propagation effects 
introduced by the nonlinear diffusion. For example, the behavior at the critical mass is different 
(PHI E] versus jU]) as well as the blow-up dynamics [23 E31 E21 E]- To treat the M. 2 case with our 
methods would require more sophisticated analysis techniques in the global existence argument (see 
Remark [9] in ^2] for a discussion) , as errors introduced by the pseudo-differential approximations 
must be controllable in terms of the Boltzmann entropy. There are also advantages in the blow- 
up argument due to the spatial localization provided by the finite speed of propagation in d > 3. 
Despite these potential difficulties, we expect analogous results to hold in M 2 . 
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1.0.1 Notation 



In what follows, we denote ||u|| p := IMIlp^) where L p (R d ) := L p is the standard Lebesgue space. 
We define the homogeneous and inhomogeneous Sobolev spaces as the closure of the Schwartz space 
under the norms \\u 

||2 := j |£|2« | fl ^|2 d £ and || u ||2^ . = + ^2y | fi ^|2 d ^ where denotes 

the Fourier transform of u(x). We will often suppress the dependencies of functions on space and/or 
time to enhance readability. The standard characteristic function for some S C R d is denoted I5 
and we denote the ball Br(xq) ■= {x G M d : \x — xq\ < i?}. In addition, we use f fdx := J Rd fdx, 
and only indicate the domain of integration where it differs from M. d . We use Af to denote the 
Newtonian potential: 

f ^log|x| d = 2 

Af(x) = < r(d/2+i) 1 12-d . ^ o 

In formulae we use the notation C(p,M, ..) to denote a generic constant which depends on pa- 
rameters p,M,..., which may be different from line to line or term to term in the same formula. 
In general, these constants will depend on more parameters than those listed, for instance those 
which are fixed by the problem, such as a(x) and the dimension, but these dependencies are usually 
suppressed. We use the notation / & . g to denote / < C(p,k,..)g where again, dependencies 
that are not relevant are sometimes suppressed. 



1.1 Background 



The local theory for (1.1) is studied in [T]. Here we simply discuss the results of that work, which 
follows closely the work of |11| l4"0"t EJ HI [2] . We begin with the definition of weak solution, which is 
stronger than the concept of distribution solutions. The main purpose of this definition is to ensure 
that weak solutions are unique. 

Definition 1 (Weak Solution). A function u(t,x) : [0,T] x M. d — > [0, 00) is a weak solution of 
( fTTj ) if u € L°°((0,T) x R d )nL°°(0,T,L l (M d )), u m e L 2 (0, T, i/ 1 (R d )), itVc G L 2 ((0,T) x R d ), 
u t 6 L^T,^- 1 ^)), and for all test functions (j) G H l {R d ) for a.e t G [0,T], 



{ut(t),<p) 



(-Vu m (t) + u(t)Vc(t)) ■ dx, 



(1.3) 



where c(t) is the strong solution to the PDE —V • (a(x)Vc(i)) + , y(x)c(t) = u(t) which vanishes at 
infinity. 



We state the following theorem summarizing the local theory of (1.1), developed in pQ as well 
as HQlElE]. 

Theorem 1 (Local Existence and Uniqueness). Let d > 3, let a(x) G C 1 be strictly positive such 
that both a{x) andVa(x) are bounded, letj(x) G L°° be non-negative anduQ G L^(M rf ; (l + |x| 2 )cfcc)n 
L°°(]R rf ). Then there exists a maximal T+(uo) > and a unique weak solution u(t) to (1.1) which 
satisfies n G C([0, T]; L l + (R d ; (1 + \x\ 2 )dx)) n L°°((0, T); L°°(R d )) for allT < T + (u ) and u(0) =u . 
Additionally, J-(uq) < 00 and we have the energy dissipation inequality, 



u(s) 



V- 



m 



u m ~ L (s)-Vc(s) 



dxds < J-"(«o). 



(1.4) 



We also have the continuation criterion: if 



lim sup \\(u(t) — k) + \\i = 0, 
k ^°° te[o,T + (u )) 



then necessarily T + (uq) = 00 and u(t,x) G L c 
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The results regarding the critical mass in the constant-coefficient case are summarized in the 
following theorem. 

Theorem 2 (Critical Mass [HIE])- Suppose a(x) = a and j(x) = 7 are both constants. Then the 
sharp critical mass satisfies, 

M _ ( 2a \ d ' 2 
Mc -\{m-l)C*c d ) ' 

and ifu(t, x) is a weak solution to ( 1.1 ) with M{uq) < M c , then u(t, x) exists globally, e.g. T + (uq) = 
00, and we have u(t,x) G L°°(M + x M d ). Conversely for all M > M c there exists a solution to 
(1.1) which blows up in finite time with M(uq) = M. Here C* is the optimal constant in the 
Hardy-Littlewood-Sobolev inequality (1.6) below and c^ is the normalization factor in the Newtonian 
potential: 

C - r ^ + 1 ) (1 5) 

° d - d(d- 2)^/2- (i ' 5j 

In [9], Blanchet et. al. exhibit a unique family of stationary solutions to the scale-invariant 
problem, based on which we will construct barriers used in the proof of finite time blow-up in 
section 3. 

Theorem 3 (Stationary Solutions to Scale- Invariant Problem [9]). There exists a non-negative, 
radially symmetric, non-increasing function V{x) supported in the ball of radius one with ||V[|i = 

(m-i)C c d J which is the unique solution (up to L 1 scaling and translation) of 

AV m = V • {VVN*V). 
Remark 1. Note that if a > and V = a d / 2 V, then 

AV 2 ~ 2/d = -V • (VVAf * V) 
a 

and in light of Theorem [2] above, V are the unique (up to L 1 scaling and translation) stationary 
solutions to the problem 

u t + V • (uVc) = Au 2 - 2 / d 
—aAc = u. 

In [9], it is also shown that these stationary solutions are the unique extremals of the following 
Hardy-Littlewood-Sobolev type inequality. 

Theorem 4 (Sharp Hardy-Littlewood-Sobolev Inequality |9J). There exists some optimal C* > 
depending only on the dimension such that for all f G L\ n L m , 

f(x)f(y) \x - y\ 2 - d dydx < C4f\\ 2 ~ m \\f\\Z (1.6) 

and equality is achieved if and only if there exists «o G M, xq G M. d , \q G (0, 00) such that 

p f \ OiQ fx - x 



We will also use the following more general inequality, although we have no need of the sharp 
constant, 

J J fix)g(y)K(x-y)dxdy< p , q> t\\f\\p\\g\\q\\K\\ L t, 00 . (1.7) 
Here, || • ||it,c» denotes the weak I) space. 
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1.2 Summary of Results 

We now state the main results. The higher regularity on a{x) and j(x) is assumed so we may use 
standard symbol classes. It is likely that significantly less regularity is required for the same results. 

Theorem 5 (Global Existence at Subcritical Mass). Let d > 3 and a(x) G C°°(R d ) be strictly 
positive and 7(x) G C°°(M. d ) be non-negative such that D a a and D a ^/ are bounded for all multi- 
indices a. Then we may estimate the critical mass as 

M c =f 2 / nf ^; (x) V /2 , (i.8) 

\{m- l)C±c d J 



and any weak solution u(t) to (1.1) with M{u) < M c exists globally and u(t) G L°°(IR + x K ). 



Remark 2. Using the methods of [2], Theorem [5] can be extended to cover more general 
filtration equation- type nonlinear diffusion on the RHS of the PDE for u(t) (e.g. AA(u), A G C 1 
and non-decreasing with < lim inf z ^oo A' (z)z 2 l d ~ l < oo, instead of simply Au 2 ~ 2 / d ). 

The next result is a significant refinement of the previous theorem, providing a strong blow- 
up criterion for solutions of arbitrary mass which requires them to concentrate a specific amount 
of mass at blow-up time. This theorem indicates that blow-up is indeed necessarily a localized 
phenomenon and ultimately does not depend on the global or average properties of a(x). Theorem 
[5] is a direct corollary of this result, however we have chosen to state and prove them separately as 
the proof of Theorem [5] is elementary once the proper machinery is in place whereas the proof of 
Theorem [6] is significantly more technical. The methods we employ to prove Theorem [6] rest on a 
careful estimate of the free energy and are based on iterating a refined version of arguments found 
in |10j . In order to do so, we employ a delicate geometric decomposition and the concentration 
compactness principle [32j [M]. The proof of [TO] is also a type of concentration compactness, 
and Theorem [6] represents the natural extension to more sophisticated concentration compactness 
methods which had not (to our knowledge) been previously used in the study of Patlak-Keller-Segel 
models. 

Theorem 6. Letd>3 and a(x) G C°°(R d ) be strictly positive and -y(x) G C°°(R d ) be non-negative 



such that D a a and D a j are bounded for all multi-indices a. Let u(t) be a solution to (1.1) which 
blows up at T+(uo) G (0,oo]. For any sequence of times tk /* T+ we may extract a subsequence (not 
relabeled) and find a sequence of points xu G M. d such that 

limsuplimsup / u(tk,x)dx> I — — — ) liminf a(xk) d ^ 2 - (1-9) 

r->0 fc->oo JB(x k ,r) V ( m — l)^* c d / k ^°° 

Remark 3. For any sequence of times tt /*" T+ we may extract a subsequence and a finite Borel 
measure n* such that uft^dx n* . If x^ is bounded then we may extract a further subsequence 



such that limfc-^oo xu = xq. In this case, (1.9) implies 

In the homogeneous case, {xk} is always bounded due to the standard estimate on the second 



moment. In the case of (1.1), we suspect it may be possible that {xk} be unbounded (i.e. mass 
will escape to infinity) if the infimum of a(x) lies at infinity. Although if T + < 00 we expect (but 
cannot confirm) that {xk} are always bounded. 
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Remark 4. For the homogeneous problem in 2D an analogous, but stronger, result on the Torus 
by Senba and Suzuki [38] shows that at blow-up time the measure n* consists only of atomic 
parts with at least critical mass and an absolutely continuous part which is smooth away from the 
concentrations (see also the e-regularity results [HIE]). However, the methods of [381 HD H2] are 
completely different from ours, which are purely variational. 

For solutions with precisely critical mass, called threshold solutions, Theorem [6] implies a very 
specific blow-up structure: in order for a threshold solution to blow up, it must concentrate all 
of the mass into a single point where the minimum of a(x) is achieved. For radially symmetric 
threshold solutions, this is enough to allow us to prove global existence and uniform boundedness 
(Theorem uy. This is in contrast to the M 2 case, where threshold solutions with finite second 
moment are global but aggregate as t — > oo |10j . This is not the case in d > 3, where critical 
solutions are uniformly bounded and cannot concentrate past a certain amount [H]. In [35], a 
symmetrization inequality from [30] (see also [H]), allows one to extend the result for radially 
symmetric solutions to all solutions. Currently no such symmetrization inequality is available in 
the inhomogeneous case, so we are restricted still to the radially symmetric result. The result is 
proved using a mass supersolution, however, it is important to note that since radial monotonicity 



is not generally conserved by (1.1), unlike the homogeneous case, a mass supersolution is insufficient 
to deduce Theorem [7j hence the necessity of also using Theorem [6] (see ^4] for further discussion) . 

Theorem 7 (Global Existence at Critical Mass). Let d > 3 and a(x) G C°°(M. d ) be strictly positive 
and 7(x) G C°°(IR rf ) be non-negative such that D a a and D a j are bounded for all multi-indices a. 
Suppose additionally that a(x) and "f(x) are radially symmetric. Then any radially symmetric weak 
solution u(t) to pTTj) with M{u) = M c exists globally and u(t) G L°°(R + x R d ). 



The following summarizes our finite time blow-up results, which demonstrates that under certain 
hypotheses, Theorems [5] and [7] are sharp. 

Theorem 8 (Finite Time Blow-Up). Let j(x) = and let a G C 1 (IR <i ) be radially symmetric, strictly 
positive and such that both a and Va are bounded. Suppose also that a(0) = mina(x) and that there 
exists a neighborhood \x\ < 5o such that a(x) is radially non- decreasing. Then for all M > M c , 
there exists a solution u(t) with M(u) = M which blows up in finite time, e.g. T + (u(0)) < oo. 

Remark 5. The requirement that a and Va be uniformly bounded are only used to satisfy the 
hypotheses of Theorem [TJ which ensures we have a well-understood local existence, uniqueness and 
stability theory. 

Naturally, blow-up solutions constructed in Theorem [8] are required to concentrate a sufficient 
amount of mass near where that minimum is achieved. Exactly how concentrated the initial data 



is required to be is characterized by (1.10) in the following proposition (which implies Theorem 
[8]), which requires at least part of the initial data be more concentrated than a particular rescaled 
extremal of the sharp HLS (Theorem [4]). Theorem [8] is proved by comparing the true solution 



against a barrier, and (1.10) below is the requirement that the solution and the barrier are ordered 
at time zero. Remark [6] clarifies how (1.10) requires uq to concentrate around where the minimum 



is achieved at least when M(uq) \ M c 

Proposition 1. Assume the hypotheses of Theorem^ hold. Let u G L\(R d ; (1 + \x\ 2 )dx) n C°(R d ) 
be radially symmetric such that M c < M(uq) and suppose that there is an Rq < 5q and an Mq with 
M c < Mo < M(uq) such that for < r < Rq we have 



|x|<r \ Ro J \RoJ \Mo/ J\ x \<r 
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Then the weak solution u(t) associated with uq blows up in finite time. Moreover, if we define 
H := M c /Mq < 1, we have the following estimate of the blow-up time: 

where a is the surface area of the unit sphere in M d . 
Remark 6. In light of Remark [TJ ( |1.10 ) implies that 



Mo(^V /2 < f u ( x )dx. 

V a (°) / J\x\<R 

Hence, in order to construct blow-up solutions with M{uq) \ M c , in general we need to choose 
increasingly concentrated initial data by sending Rq — > 0. 

Remark 7. Arguably, the most natural inhomogeneity in the PDE for u(t, x) is a mobility (which 
weights the cost of mass transport) that depends on space, that is b(x) > such that 

ut = V ■ {b{x)uV^u m - x - Vc) 

-V • (a(s)Vc) + 7(k)c = u. K ' ' 



In this case the free energy is still given by ( |1.2[ ); the gradient flow structure is only changed via 
the metric. The Wasserstein metric would be replaced by a generalized version (still based on 
L 2 ) that accounts for the inhomogeneity in the cost of mass transport. As Theorems [5] and [6] are 



proved using essentially only (1.2), the results still hold for (1.12) without any relevant change in 
statement or proof. The result of Theorem [8] will also hold under suitable monotonicity hypotheses 
on b(x) and the proof is only a minor modification of the current one. Other possible variants of 



(1.1) could involve mobilities which are different in each term, however this would probably change 
the variational structure and the methods we employ to prove Theorems [5] and [6] would no longer 
apply. Regardless, this is a very different model from a physical standpoint, as this corresponds to 
adjusting the organism response to stimulus dependent on space. 

As mentioned above, Theorem[8]and Proposition [T] also provide new results for the homogeneous 
problem, a(x) = 1, j(x) = 0. The virial method used in \9] proves blow-up for all solutions with 
negative free energy, without the need for radial symmetry. On the other hand, Theorem [8] and 
Proposition [T] require radial symmetry but do not require any assumptions on the free energy. 
Indeed we have the following corollary of Proposition [TJ which in particular, shows the existence of 
solutions with arbitrarily large initial free energy that blow up in finite time. 

Theorem 9. Let uo £ L\(R d ; (1 + \x\ 2 )dx) D C°(R d ) be radially symmetric with uq(0) > and 
satisfy 

d/2 

M(u ) > ' 



(m - l)C*c d/ 

Then the weak solution u(t) associated with uq of the scale-invariant problem (a(x) = 1, 7(x) = 0) 
blows up infinite time. In particular, for every Fq > 0, there exists a solution u(t) with F(u(Q)) > Fq 
which blows up in finite time. 

Remark 8. It would be interesting to determine if there exist solutions which have positive energy 
at blow-up time. Such a result would require deducing more precise information about the blow-up 
structure. 
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2 Global Existence 

The proof of Theorems [5] and [6] hinge primarily on providing a precise decomposition of the potential 
energy, 

u(x)c(x)dx, 



into a leading order critical part and another part that is subcritical. The use of pseudo-differential 
operators for this purpose is the subject of the following section. 

2.1 Approximate Inverse of Chemo-attractant PDE 



We use the standard symbol classes studied in for example |39j . summarized in the following defi- 
nition. 

Definition 2 (Symbol Class S s , s E M). Suppose € C°°(IR^ x Rg) satisfies 

^Kx,o\<fi, a (i + ieir H , 

for multi-indices a, f3. Then we say both b and the associated pseudo-differential operator (^DO) 
Tfe defined by 

are in the symbol class S s and say the symbol b or the operator Tb are of order s. We also denote 
6(a?,£) = sym(T fe ). 

Notice that with this definition S Sl C S S2 whenever s\ < S2- Also, since the symbols are required 
to be smooth, these operators do not carry too much low-frequency information, unlike multiplier 
or symbol classes that allow singularities at the origin. For the standard relevant facts regarding 
these symbol classes, such as the symbolic calculus, localization estimates, boundedness on Sobolev 
spaces and singular integral representations, see Chapter 6 of [39] . 

Consider the PDE 

Lc := -a(x)Ac - Va(x) ■ Vc + 7(z)c = u. (2.1) 
By definition, L is a pseudo-differential operator in S 2 : 

Lc = j^j- 2 J ( 7 (x) + a(x) |e| 2 - i£ • Va(x)) 
Consider the approximate inverse of L, the S~ 2 class ^DO 

a u-= 1 ! mme^ 

H ' (2vr) d / 2 7 a{x)\i\ 2 -ii-Va{x) + 1 {x) 

where $(£) = E[j=i ^(C?) 0(0 a smooth function such that < < 1 which is identically 
one for \t\ > 1 and vanishes in a neighborhood of zero. We remark that if j(x) is strictly positive, 
we do not need the cut-off $(£). By the symbolic calculus [Chapter 6, Theorem 2 |39|], 

A H Lc = c + T E c, 
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where the operator Tg £ S 1 and the associated symbol E has the following asymptotic expansion 



£-$(0-1+ £ (2 " } , '"' ^yml.W/i^sy.nd). 

|o|>l 



in the sense that the error in truncating the series for N > |a| is a symbol of class S~ 



-1-JV 



Although both Ah and Te are bounded operators from LP to itself for 1 < p < oo (Chapter 
6 [39]), if 7(x) is not strictly positive then -L -1 , the true inverse, is not. Necessarily, the low- 
frequency portion of is still present implicitly in Tec. Instead of being bounded on LP to itself, 
satisfies the following: for 1 < q < d/2, d/(d — 2) < p < oo and | + ^ = ^ 



l L 1u IIp ^ IN 



(2.2) 



This can be seen at least formally by multiplying both sides of (2.1) by c^ d 2 ) p / d 1 integrating, 



and applying the homogeneous Sobolev embedding (in the homogeneous case, (2.2) is the Hardy- 
Littlewood-Sobolev inequality). 

We may formally write down the operator Ah as a singular integral operator by interchanging 
the integrals: 



A H u(x) 



1 



(2tt) c 
1 

(2tt) c 



u{y) 



$(£)e 



o(x) |£| - if • Va(x) + 7(z) 



<%df 



$(£)e 



i((x-y) 



a(x) |f| J - if • Vo(i) +7(z) 



u(y)K H (x,y)dy. 



The integral for Kn(x,y) is not absolutely convergent so we cannot naively apply Fubini's theorem 
in the above computation rigorously, but it can be justified by a standard limiting procedure, as 
in |39j . The key technical lemma for the proof of Theorem [5] is the following characterization of 
Kn{x,y). 

Lemma 1 (Asymptotic Expansion for Kff(x,y)). Let Kff(x,y) be defined as above by the condi- 
tionally convergent integral 



Ku{x,y) 



1 



(2vr) d 7 a(x)\t\ 2 - if • Va(x) + 7 (a 
Then we then have the following asymptotic expansion which holds uniformly in x G 

KH{x,y) = -^-r \x - y\ 2 ~ d + o(\x - y\ 2 ~ d ) as y -)■ x, (2.3) 
a{x) 



with Cd given above by (1.5). Moreover, recall that for all 5 > and N > (see for example pg 235 



\k h {x,: 



< 



s,n \x - y\ 



-N 



\x — y\ > 5. 



(2.4) 
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Proof. The bound (2.4) is a standard consequence of Ah 6 S~ 2 . Such localization should not be 
surprising since the the low frequency contribution of L~ l is not included in Ah due to the cut-off 
Hence, we focus on (2.3). Note the trick 



1 

D 



e- tD dt. 



Hence, 



Kh(x,v) 



(2tt) c 



^^ e i (i x -y+ t ^ a ( x ))- ta ( x M\ 2 - t 'r( x ) ( ] / ^ ( ] i f i 

) poo 

e -ty{x) / 0(^.) e <O(*i-»i+*«« i *(*))-to(*)^ d ^ dt . 



Now define the complex change of variable Zj 



\ia\X)) 5 3 Z 2 (ta(x))V2 > 



4ta(x) 



(ta) 1 / 2 



R-G Zj \ _ - - 

(ta(x))V2 ) e ^ 



4ia(x) 



(ta(x))V2 



where Tj is the contour |lmzj 



Xj-yj+td x ,a(x) ' 



-td x .a(x) "I 
2(to(x)) T 72 /• 

Applying the above change of variables to the expression for Kn(x,y) implies 



Kh(x,v) 



-(x—y)Va{x)/(2a(x)) poo 



t 



-d/2 



{2ir) d {a{x)) d / 2 J 
We make the following additional change of variables, 



-t'y(x)-t ' Va{x >' |j ~ y| , 

1 'W 1 4a(x) 4ta(i) 



t 



which then yields 



c 



\x - y\ 



(4to(x))V2 

In terms of £, Kn(x,y) can be now written as 



e -(x-y)Va(z)/(2a(x)) 

W(2a(x)) |X_2/r 



4a(x)C 2 ' 



and eft 



F ~ j/l 
'2a(x)C 



-d(. 



-d-2, 



3=1 



x — y\ { \Wa(x)\- . , i 



Due to the smoothness of and a and the strict lower bound on a we have the uniform (in x) 
convergence of the integral (note that here £ is fixed and Zj is the complex integration variable) 



lim/ i (t(C)) = lim / 



\x — y\ 1 



e -*idZj = vr 1 / 2 . 
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Similarly we also have 



lim 



-d-3 



g-»r I |Va(g)l 



<Z/2 



uniformly in i £ ff 1 due to the uniform continuity of Va(x),a(x), and 7(x) as well as the strict 
positivity of a. Recalling elementary facts about the Gamma function, 



OO 1 



d 



Hence, 



r(l + 1) 



d(d - 2)vr d / 2 a(x) 



2 /d 
d(d - 2) V2 + 



I |2-d , / i i2— d 

\x-y\ +Oy^ x [\x-y\ 



and (2.3) is proved. 



□ 



The error term Tec in the approximate inverse can be controlled by the following lemma. Here 
we take advantage of the smoothing nature of T E £ to show that in the potential energy, this 
term is subcritical in the sense that the effective power of || "ix 1 1 associated with the term is strictly 
less than m. Naturally, one must eventually use (2.2) in order to prove this lemma. 

Lemma 2 (The error term is subcritical). Let d > 3, suppose u G L 1 (M. d ) D L m (M. d ) and c is the 



strong solution of (2.1) which vanishes at infinity. Then, 

J uT E cdx < \\u\\j 



2-20H \\26 
l"llm > 



(2.5) 



for some 0<29<m = 2- 2/d. 
Proof. Define 



4d 



2d - 3 



Define the standard Fourier multiplier (V) s f := ( 1 + |£| J /(C)- Since the multiplier (V) 1 ' 2 is 
self- adjoint we have, 



uT E cdx 



< ikv)- 1 / 2 ^ \\{V) x l 2 T E c 



p-i 



Define 



1 1 2 2d + 5 
q p d 4a 



Definition ^implies that since T E E S' 1 , we also have T E G S' 1 / 2 . Hence, T E is a bounded operator 
L p to W 1 ' 2 ' p [Chapter 6, Proposition 5 |39j]. Using this and the elliptic LP estimate (2.2) we have 



One may easily verify that 



\m l l 2 T E c\\ v <\\c\\ v <\\u\\ q . 



1 _ p- 1 1 

<7 p 2d' 
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and therefore the inhomogeneous Sobolev embedding theorem implies 

'1/2, 



[|<v>-^«|u <im 
p-i 



Hence, we see the relevance of q as we have in total, 

uTeccIx 



< llu" 2 



In order to interpolate between L 1 and L m , we need q < m = 2 — 2/d, which follows easily from 
d > 3. Then, for 

(2d-5)(2d-2) 



Ad(d - 2) 



G(0,1) 



we have, 



uTEcdx 



< Ilu||^||ull2?. 



To prove subcriticality it remains to confirm that we have 29 < m, which again follows from 
d> 3. □ 

2.2 Proof of Theorem g] 

In this section we complete the proof of Theorem [5] 

Proof. We prove Theorem [5] by producing a uniform in time bound on the entropy (which is basically 
just the L m norm). This in turn proves that the solution is uniformly equi-integrable and Theorem 
[TJ completes the proof. 



By the energy dissipation inequality (1.4) and the definition of Ah, 

/u m dx / uAnudx H — / uT^cdx < T{uq). 
2 J 2 J 

By (2.5) (Lemma [2]) we then have, for some 0<29<m = 2 — 2/d, 

^-j J u m dx - i J uA H udx - CWuf^Wuf^ < T(u ) 



m 



Using Lemma [TJ for every e > 0, we may choose a 5 > such that 

( _Cd_ \ u(x)u(y) 

*\x-y\<5 



uAuudx < 



a(x) 



x 



y\ 



id— 2 



dxdy + C(6)\\u\\i, 



where Cd, given by (1.5), is the normalization constant in the Newtonian potential. Hence, by the 



sharp Hardy-Littlewood-Sobolev inequality (1.6) we have, 

C* ( c 
2 \mma(x) 



1 



m 



1 



+ e 



1 2/d 



<C||n||?- 2e h||^ + C(5)|| u ||^+^ ). 



1 29 



If M(u) < M c with M c given by (1.8) then we may choose e sufficiently small such that the first 

t . This 

□ 



term is positive. Since m > 29 this then implies a global uniform- in-time bound on \\u 
in turn implies a global L°° bound on u(t) by the continuation criterion in Theorem [TJ 



Remark 9. It appears the proof of Theorem [5] would require some non-trivial refinement in order 
to treat the M 2 case. Naturally, the asymptotic expansion of Kfj(x,y) along the diagonal x ~ y 
would need to be refined in order to capture the logarithmic singularity accurately (Lemma [TJ). A 
likely more delicate task would be adjusting Lemma [2] to yield an error which is subcritical relative 
to the positive part of the entropy f u(\ogu) + dx. 



14 



2.3 Proof of Theorem M 

Like Theorem [5j the proof of Theorem [6] is essentially a precise estimate of the potential energy 
which allows the free energy to uniformly control the L m norm of the solution. In Theorem [5j 



after the results of 6 2.1, the assumption of subcritical total mass was used to make this estimate. 
Theorem [6] refines this result to show that the free energy uniformly controls the L m norm unless 
a critical amount of mass concentrates into a single point or at least along a subsequence of points 
escaping to infinity. The main idea of the proof is to combine the related result for critical mass 
solutions in 2D found in [10J with a more intricate concentration compactness-style decomposition. 
It is a little technical, but as with many concentration compactness arguments, the key idea is 
simple: if the solution does not concentrate a critical amount of mass at blow-up, then we may 
divide 3R d into a finite number of domains, each containing subcritical mass at blow-up. In each 
region we may apply reasoning similar to Theorem [5] to control the corresponding contribution to 
the free energy, noting that long-range interactions between different regions are "sub-critical" since 



the leading order term in the approximation of the potential energy deduced in {2.1 will not be 
present. 

Proof. We suppose Theorem [6] is false and go towards a contradiction. For the duration of the 
proof, the notation B(x,r) will denote the closed ball with center x and radius r. Moreover, let 
a = min a, A = max a and define 



M* 



[m - l)c d C+ 



d/2 



The assumption that Theorem [6] is false is equivalent to the existence of a time sequence t k — > T + 
and e > such that for all sequences {xk}^ =1 C M. d , 

lim sup limsup / u(t k ,x)dx - a(x k ) d / 2 M* < -2e. (2.6) 

r-»Q fc^oo \JB(x k ,r) J 

Let u k {x) := Uk(tk,x). The purpose of the additional factor of two is purely cosmetic and will be 
apparent in what follows. Since the proof mostly rests on an estimate of the potential energy, define 
the following (since L~ l is self-adjoint) 



PE(f,g) := l -j fL-'gdx = l - f gL~ l fdx. 



The proof begins with the following equivalent version of the concentration compactness lemma 
found in [32] 

Lemma 3. Let be a sequence of non-negative functions in L l (M. d ) with = M S (0,oo). 

Then for all 5 > 0, there exists a subsequence (not relabeled) of {u^} such that the following holds: 
there exists an L < oo such that there exists sequences x l k C ~R d ,R l G (0, oo) and u l k ,ek,u[ 6 L 1 
(which all depend on 5) for all I G {1, L} which satisfy 

L 

u k = ^2 u i + u k + e k, (2-7) 
i=i 

with the following properties: 

(i) The supports of u l k , e k and are all disjoint. 
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(ii) suppu^, C B(x l k , R l ) and lim^oo \x k — x l A = oo if I ^ m. 
(Hi) limsup fc _ > . 00 ||e fc ||i < 5. 

(iv) For all I, lim^^ dist (suppu^suppuj^) = oo. Moreover, uY is vanishing in the sense that 
for all R > 0, 

lim sup / u k (y)dy = 0. 
k ^°° x m d Jb(x,r) 



Proof. The proof proceeds by iterating Lemma II. 1 in [32]. As the starting point of the algorithm, 
define p\ = u&. Now consider the inductive step as follows. Suppose that u l k which satisfy the 
conclusions of Lemma [3] have been determined for < / < j and suppose that v l k are a sequence 
which satisfies the properties: 

(a) limsup^oo < 2- [ - 1 S 

(b) For all < I < j, the support of vi are disjoint from uV? for all m and v™ for all m 7^ I. 

We are also given p k := u k — Ylo<i<j n L + v i- By Lemma II. 1 in [32] one of three possibilities occurs 
up to extraction of a subsequence (not relabeled) . If pi is tight up to translation then there exists 
a sequence x k and a R 3 such that 



/ 



p{(y)dy < 2-35. 

\B(x{,W) 



Hence, define u{ := and e k := Pk\ R d\ B ^i ^ + J2o<Kj v l- We ma y terminate the 

algorithm here and set L = j. If instead p k is vanishing then simply define u k := p k and = 
^2o<l<j v k- Terminate the algorithm here and set L = j — 1. In the case of dichotomy, decompose 
pi as in Lemma II. 1 of [32] as 

A = pk 1 + pf + v l 

with limsup^^ \\v k \\i < 2~ J ~ 1 5. The following additional properties are satisfied: there exists a 
sequence xj, and an R 3 > such that pfe 1 = p{.\ B ( x 3 and a sequence R k > with lim^oo R k = 00 

such that f? k = Pk\^d\ B ( x j Rk y This latter fact is what ensures (ii) and the first part of (iv). Take 
u\ = pfc 1 (with the associated x k and R J ). If limsup^^ ||Pfc 2 ||i < 5/2 then terminate the algorithm 
here with L = j and define = p* k 2 + J2o<i<j v k- ^ no ^ continue the algorithm with p^ +1 = pj; 2 . 
This termination condition ensures that the algorithm terminates in at most 0(5~ 1 ) steps. □ 

We will fix 5 = 5(a, e) later depending only on e, a and the total mass M. We apply Lemma 
[3] to Uk := u(tk), tk /* T + with this 5. This supplies the first step in the decomposition, which 
must be refined further. The following geometric decomposition lemma refines the decomposition 
by subdividing the u l k into pieces of subcritical mass. The authors would like to acknowledge Jonas 
Azzam for assisting in the proof, which is based on an iterative application of the Besicovitch 
covering theorem |35j . We postpone the proof until after the proof of Theorem |6j 

Lemma 4. Let Uk = Uk(tk) as given above. Then for any 5 > there exists a collection of 
N = N(5) < 00 disjoint closed balls with centers x\ C M. d and radii rj > 0, 1 < j < N which satisfy 
the following properties for all k sufficiently large, 
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/ \d/2 

(1) For allj, J B ^ r ^u k (y)dy < ^liminf fc _ ) . 0O mmy_^ | <r . a (a?) J M* - e. 



"7 »»' 



There exists a constant cq > independent of k such that mini<j 3 -<jv 

We proceed to prove Theorem Theorem [6j with the decomposition obtained in Lemma [4} Let us 
denote B k j = B(x 3 k , Vj) and E k := W i \ (uj^-Bfcj). These balls allow us to decompose the potential 
energy as follows, 

PE(u k ,U k ) = ^ PE ( u k\B k j,U k \ Bk J + ^2 PE ( U k\B[^ U k\E k ) + PE(u k \ Ek ,U k \ Ek )- (2.8) 

l<j,i<N l<i<7V 

First consider the last term, using the pseudo-differential approximation, 

PE(u k \ Ek ,u k \ Ek ) = \ u k (x)u k (y)K H (x,y)dxdy - u k (x)T E L~ l (u k \ Ek ){x)dx. (2.9) 

J JE k xE k JE k 

The error term is controlled by Lemma [2j 

/ u k {x)T E L- 1 (u k \ Ek ){x)dx < C E \\u k \ Ek \\™ < C E \\u k \\™, 
JE k 

where C E is the implicit constant in (2.5) times 1 1 1 1 1 261 - Now we turn to the leading order term in 
(2.9). Using Lemma [TJ for all e > 0, there is an n > such that 



u k (x)u k (y)K H (x,y)dxdy = / u k (x)u k (y)K H (x,y)dxdy 

E k xE k J JE k xE k n{\x-y\<ri} 



+ u k (x)u k (y)K H (x,y)dxdy 

J JE k xE k r\{\x-y\>ri} 

<{c d + e) / — -u k {x)u k (y) \x - y\ 2 ' d dxdy 

J JE k xE k n{\x-y\< V } a{x) 

+ sup sup \Kn(x,y)\ 

x,y£M. d \x— y\>ri 

Denote 

K(e):= sup sup \K H (x,y)\ ||n||f. (2.10) 

x,y£R d \x—y\>r/ 

For this step it suffices to take e = 1. Focusing now on the first term, cover M rf by cubes of unit 
volume, denoted {Qi}^i = Q- By construction (Lemma [4] (2) and Lemma [3] (iv)), there is a ko 
sufficiently large such that 

sup sup / u k \ Ek dy < 45. (2-11) 

k>k Q&QJQ 



Hence, by (|1.7|) and Holder's inequality, 

-u k (x)u k {y) \x - y\ 2 ~ d dxdy 



E k xE k 



n{|x-y|<??} a \ x ) 



V) / / — ( -x u k (x)u k (y) \x - y\ 2 d dxdy 

' 'QiXQjnEkXEknUx-y^rj} a { x ) 



<r~ \^ ll l | ■ 1 /rf i ■ i iiV^ii i nm/2n i nm/2 

^ 2^ iPfelEfcnQilli iPfcl^nQjIli ||«Jfe|B fc nOill m ' ||«fc|£ fe nQ 3 - || m 



Qi,Qj£Q 
< K x 5 2 / d \\u k " m 



rn > 
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for some constant K\ > which depends only on dimension and a. Putting the previous estimates 
together we have shown that, 

PE(u k \ Ek ,u k \ Ek ) < K 1 5 2 / d \\u k \C + C E \\u k C + K(l). (2.12) 

Let us now handle the second term in ( |2.8[ ). Using the pseudo-differential approximation, 

PE(u k \ Uj B kj ,u k \E k ) = / u k (x)u k (y)K H {x,y)dxdy - I u k {x)T E L' l (u k \ Ek ){x)dx. (2.13) 

The error term is controlled by Lemma [2j 

/ u k (x)T E L- 1 (u k \ Ek )(x)dx < C E \\u k \ UjBk J e m \\u k \ Ek \\ 6 m 

JVJjB k j 



< C E \\u k 



{26 



The leading order term in (2.13) is estimated using Lemma [TJ which implies for all e > 0, there is 
r] > such that 



u k (x)u k (y)K H (x, y)dxdy 



u k (x)u k (y)K H (x, y)dxdy 

UjB k ^xE k r\{\x-y\<ri} 

+ I I u k (x)u k (y)K H (x,y)dxdy 

IUjB k}:j xE k n{\x-y\>Ti} 



< (c d + e) 



yjjB k ^xE k r\{\x-y\<r,} a { x ) 



1 -u k {x)u k (y) \x - y\ 2 d dxdy 



where K(e) is defined in (2.10). For this step it again suffices to take e = 1. By (1.7) and Holder's 
inequality, 



-^-u k (x)u k (y) \x - y\ 2 d dxdy < \\u k \ Ek \t\ d \ 
UjB kd xE k n{\x-y\<r)} a \ x ) 



u k \\Z- 



Hence, the estimates put together imply, 



N 



^ PE(u k \ BhJ , u k \ Ek ) < K 2 5 l / d \\u k \\2 + C E \\u k C + K(l) 



(2.14) 



for some fixed positive constant K 2 = K 2 (d, M, a) 



Now we turn to the first term in (2.8). Firstly, if j ^ i then by arguments similar to above, 
using that the balls B l k ■ are disjoint for k sufficiently large, 



PE{u k \ Bk v u k \ Bkl ) < C E \\u k \\ 2 ^ + K 3 



(2.15) 



for some K% > which depends on M, d, a and the minimal distance between B k j and B k ^ which 
is bounded below by Lemma [4] (3). Hence, consider the case i = j and for notational simplicity 
refer to B k := B k) j. Estimate now as above, 

PE(u k \ Bk ,u k \ Bk ) = / u k (x)u k (y)K H (x,y)dxdy - u k (x)T E L~ 1 (u k \ Bk )dx. 

J JB k xB k JB k 
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By Lemma [2j 

/ u^x^eL' 1 (u k \ Bk )dx < C E \\u k \ 

JB k 

By Lemma [TJ for all e > there exists an r\ = 77(e) > such that 



29 
rn ■ 



Uk(x)u k (y)K H (x, y)dxdy 



+ 



Uk(x)u k (y)K H (x, y)dxdy 

(B k xB k )n{\x-y\< V } 



n k (x)n k (y)K H (x,y)dxdy 

(B k xB k )n{\x-y\> v } 



< (c d + e) 
+ K(e). 



(B k xB k )n{\x-y\< v } a ( x ) 



u k (x)u k (y) \x - y\ 2 d dxdy 



Write 



a,j := liminf min a(x). 



k— >oo x£B k 



By the sharp HLS, Theorem |4j 



Hence, 



-^—Uk(x)uk(y) \x - y\ 2 d dxdy < ^ ||u fc | B J 2/d ||-u fc |£>J™. 

(B k xB k )n{\x-y\< V } a {%) CLj 



PE(u k \ Bk ,u k \ Bk ) < ^-{c d + i)\\u k \ Bk \\ 2 i d \\uk\B k \C + C E \\uk\\m+K(e). 



(2.16) 



By Lemma[4](l) we may choose e sufficiently small depending on e and k sufficiently large to ensure 
that, for all j, 



mm 



1 



{c d +~e)\Wk\B k Al' d ) >c(a,e)>0 



i<i<iV \m — 1 2a j 

This in turn fixes K4 := K{e). Notice that depends only on dimension, a and e. 
Applying (pUb, pll, pT5l), (f2TT6l) and pTfl) to M we have, 



(2.17) 



F( u o) > F(uk) = — — 7 [ ufdx - PE(u k ,u k ) 
m — lj 



i/d 



\ u k\\m 



(N 2 + N + 2)C E \\u k \\™ - 2K(1) - K3N' 2 - NK 4 . 



Since depends only on e and K\,K% are fixed constants which do not depend on the decom- 
position, 5 = 5(a, e) can be chosen a priori sufficiently small such that there exists a constant 
K* > 0, 



T(uq) > K*\\u k \\Z -(N 2 + N + 2)C E \\u k \\Z ~ 2^(1) - N 2 K 3 - NK 4 . 



(2.18) 



Note that the choice of 5 then fixes the decomposition given by Lemma [4] (and hence N). As 
T^uq) is a fixed finite number and 29 < m, this estimate implies that [|«jb||m is uniformly bounded 
for k sufficiently large. However, note that (2.18) is generally vacuous unless ||ufc|| m is extremely 
large. Regardless, (2.18) implies that hmsup^^ ||ufc|| m < 00 which, by the continuation criterion 
in Theorem [TJ implies that u cannot blow-up at T+. □ 
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We now prove Lemma |4j 

Proof. Let u l k ,ek,u^ as given in Lemma [3] with 1 < I < L, For simplicity we only treat the 
case L = 1, as we may simply append resulting decompositions in order to treat L > 1. Denote 
h\ := u\(x + x\) and up to extraction of a subsequence, there exists some measure n* such that 
n\ ^* n* . Define the set E\ := suppra* C ^(O,/? 1 ). By the assumption (2.6), for every x € E\ 
there exists cin t x > which satisfies: 

(a) n*(B(x, r x )) < (min^i^ a(y - x l k )) d 2 M* - e, 

(b) is a point of continuity of the map r >— > n*(B(x,r)) (which is continuous a.e. since it is 
non-decreasing) . 



Denote this covering of E\ by B 1 . Notice that in order to satisfy (a) given the assumption ( |2.6| ), 
the radius of the balls must be chosen small enough such that a(x) is sufficiently close to a(y) 
for \x — y\ < r x (which requires the assumption of uniform modulus of continuity of a(x)). It now 
follows by the Besicovitch covering theorem (Theorem 2.7 [35J) that there exists a finite or countable 
set of balls B 1 C B 1 which satisfies: 



(i) l^i < Y^BeB 1 — w here P(d) is a dimensional constant. 

?i 



(ii) There exists Q{d) families of disjoint balls B\, ...,Bq,* C B 1 such that E\ C Here 



again, Q(d) is a dimensional constant. 
Since the collections Bj define a covering of E\ we have also 

Q(d) 



n*{E x ) < E E n ^ B ^ 
»=i BeBl 

which implies that one of the collections, without loss of generality suppose B\, satisfies 

® <£'">■ 

We may truncate this sum to find a finite collection of B € B\, labeled B\ i for i £ {1, Ni} which 
instead satisfies 

Now we describe how to carry the m-th step of the algorithm to the m + 1-st and justify why the 
algorithm generates the decomposition stated in Lemma [4] in finitely many steps. Suppose we are 
given a set E m+ \ C B{Q,R l ) which is defined by 



E m+ i := Ei \ Uf =1 U^L\ B\^ 



where we assume that the balls ^i?*^, 1 < i < N}~, 1 < k < m are disjoint and chosen such that 



"*< e -«) £ I 1 -2qmJ "*< Ei) = i 1 - Wi) 1 "* (E '" ) ' (2 ' 19) 
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From (2.19) it is clear that if m is sufficiently large depending only on n*(E{), then n*(E m+ i) < 5 
and, in particular, 

n*(R d \ Uf =1 U& = n*{E m+1 ) < 5. 

Hence, the primary remaining step is to see that we really can progress from step mtom+1. To 
this end, define a covering B m+1 of E m+ i by balls B(x, r x ) for all x E E m+ i which satisfy: 

(a) n*(B(x, r x )) < (min |y _ x | <riI1 a(y + x[)) d/2 M* - e. 

(b) r x is a point of continuity of the map r i— > n*(B(x,r)). 

(c) B(x, r*) n U£L X ugi = 0. 

To achieve the condition (c), generally r x must be taken small, however it is possible since the balls 
in question are closed. Given the covering B m+ , the Besicovitch covering theorem implies that 
we may choose a finite or countable subset B m+l C B m+1 which satisfies properties analogous to 
(i) and (ii) above, and moreover satisfies the additional property that B m+1 is disjoint from B\^ 
Given the sets chosen in (ii), we may again assume without loss of generality that, 

and furthermore that we may truncate this sum to be a finite number of balls labeled B™^~ with 
1 < % < N m+ i for some finite N m+ \ such that 

*(w \ Nm + 1 
2Q(d) ~ ^ n ^ ] - 

Appending j-B™" 4 " 1 } to the existing list and defining E m+ 2 by 

E m+2 := E X \ U& Bf f<> 

we see that indeed we have advanced the algorithm to step m + 1. By the argument above we saw 
that the algorithm terminates at some finite M < oo. Hence define N = Yui=i Ni and re-index 
the balls B\ i to Bj, 1 < j < N. Each ball Bj has center Xj G ^(O,^ 1 ) with radius rj > 0. We 
set x 3 , := Xj + x\. Since each rj is chosen as a point of continuity, (1) and (2) in the statement 
of the Lemma are satisfied by ut for k sufficiently large. This can be verified by continuously 
approximating 1 B ^ X r ^ and using the fact that n\ n* . and limsup^^^ ||efe||i < 5. Finally, 
observe that (3) in the statement of the Lemma is satisfied by construction and Lemma [3| □ 

3 Finite Time Blow-Up 

As discussed in the introduction, the inability to make obvious use of a virial method motivates our 
use of a barrier method based on maximum principle-type arguments. We begin with the following 
rescaling: let M c < Mo < M(uq) be as in the statement of Proposition [TJ define 

H := M c /M < 1 (3.1) 

and let 

p(t,x) := /J,u(fj}' 2/d t,x). (3.2) 
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Then M(p) = M(u)p> M c and p solves 

Pt + V • {pp x - 2 l d Vc) = Ap 2 - 2 / d 



V • (a(x)Vc) = it 1 p. (3 ' 3) 
As mentioned above, we will use mass comparison arguments involving suitably chosen barriers 



(sub-solutions) to force the finite time concentration of mass in (3.3). We define the following 
comparison function u = u(t,x) with Rq as in Proposition [TJ 

where we take R(t) as a solution to the initial value problem 

f R(t) _ M c (l-M- 2/d ) 

< ^ W ~~ a^o^iW- 1 (3.5) 
I fl(0) = Rq. 

Here a denotes the surface area of the unit sphere. Notice that R < and that R(T*) = with 

T.: = r.(Mo. Ro . d ) ^ aWM ^_. 2/ ,_ 1) <oc. (3.6) 

We define the mass distributions 

M(t,r)= / p(t,x)dx, M(t,r)= / u(t,x)dx. 

J\x\<r J \x\<r 



Notice that (1.10) is equivalent to M(r, 0) < M(r, 0), which means that the rescaled initial data is 
initially more concentrated than the barrier u on the neighborhood r < Rq. It is also important to 
note that the total mass of the barrier u is generally more than the critical mass M c but less than 
or equal to the total mass of p which itself has less mass than the true solution u. 

Suppose that M(0, r) < M(0, r) for all r > 0. We will show that this ordering is preserved up 
to the blow-up time of u or u, e.g. 

M(t, r) < M(t, r) for < t < min (T^, T + (u )) . 

See Figure 2 for a graphical depiction of the relationship between u(t, x) and u(t, x), which is very 
different from the way M(t, r) and M(t, r) are related. 



As alluded to above, in the language of maximum principle-type arguments, u plays the role of 
a subsolution in terms of the mass concentration. As u concentrates into a delta mass at t = T+, 
we must have T+(uq) < T*, which will conclude the proof of Theore m [8} The intuition for why the 
proof ultimately works is based on the fact that the rescaled system (3.3) pulls mass into the origin 
faster than the PDE that a d ^ 2 (RQ)V solves (see Remark [l]). That is, the rescaling (3.2) transfers 
the property of having supercritical mass into surplus attractive power when compared against 
stationary solutions of roughly comparable mass. The surplus attractive power is what gives us the 
ability to choose R(t) at the rate given in (3.5) and hence prove that the solution is concentrating 
fast enough to be squeezed into blow-up by the self-similar barrier u. 
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Figure 1: The mass subsolution u(t,x) and the real solution u(t,x). 

We use mass comparison arguments influenced by those found in [30] to prove that the mass 
ordering ^ is preserved (see also [5j E7J CE] ) • As mentioned in the introduction, the novelty in our 
argument lies in the nature of our barriers which can handle near-critical mass accurately up till 
blow-up time. We should mention that, due to the inhomogeneity, the mass comparison principle 
that normally holds (in the sense of aforementioned references) does not apply. 

A nontrivial complication arises at the free boundary of the positivity set of u(t). Here, classical 
regularity breaks down and the mass comparison arguments no longer provide a rigorous argument. 
To deal with this technical issue we lift to strictly positive solutions, which are smooth due to 
uniform parabolicity on compact sets. Comparison with vanishing error is proved against these 
solutions, for which the formal arguments are rigorous. Passing to the limit requires some standard 



approximation arguments regarding the stability of (3.3). In order to prove positive solutions remain 
positive, we adopt viscosity solution-type arguments similar to those used in [29] (see Appendix for 
more information). 

3.1 Preliminaries: Approximation and Regularity 



In this section we detail the important approximation and regularity properties of (1.1). These 
results are more or less expected, but since they are of independent interest and important for 
making our arguments rigorous, we include brief sketches of the proofs. 

Lemma 5. Suppose uq G C°(M d ) n L+(M. d ; (l-\-\x\ 2 )dx) and let u(t) be the associated weak solution 
which satisfies u(0) = Uq. Then for all e > 0, 

(a) Vc is continuous and bounded, and Ac is bounded on t € [0, T+(uo) — e)- 

(b) u(t,x) G C°([0,r + (u ) -e) xR d ). 

Proof. (a) The proof follows from standard elliptic regularity estimates and is omitted for brevity. 

(b) By (a), our PDE ut — Au 2 ~ 2 / rf + V • (uVc) = can be viewed as a degenerate diffusion with 
a priori given drift term c(x), which satisfies the assumptions (A1)-(A3) of DiBenedetto |21| . 
In addition, from the proof of local existence theory (see [2]), the approximation assumption 
(A4) is satisfied (i.e. u can be approximated locally uniformly with the smooth solutions u n 
with strictly positive initial data). Therefore Theorem 1.1 of DiBenedetto [21] yields (b) (also 
see Theorem 3.1 of |30j). 

□ 
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Lemma 6 (Regularity of Strictly Positive Solutions). Let u € L\(R d ; (1 + \x\ 2 )dx) n L°° (R d ) be 
strictly positive a.e.. Then for all e > 0, the associated weak solution u(t) with u(0) = uq is smooth 
and strictly positive on (0, T + (uq) — e) x R d . 

Proof. Due to Lemma [9] in the appendix, u stays strictly positive for all < t < T+(uq). Con- 
sequently u is a solution of a uniformly parabolic quasilinear PDE of divergence form, and the 
regularity of u follows from classical regularity theory [31] . □ 

Lemma 7 (Stability of the blow-up time). Let {u%} C L\(R d ; (1 + \x\ 2 )dx) n L°°(R d ) and u n (t) 



denote the associated weak solutions of (1.1) with u n (0) = Uq defined on the intervals [0,T+(mq)). 
Suppose further that 

(a) sup n (||«ol|oo + Kill) < 00 an d 

(b) u% -)■ u strongly in L 1 for some u € L\(R d ; (1 + \x\ 2 )dx) n L°°(]R d ). 



Let it(i) be the solution to (1.1) u(0) = uq; an <^ define Tq > suc/i £/&a£, 



To := sup < T G (0,oo) : liminf sup ||ti n (i)||oo < 00 



T/ien /or aZZ T < Tq, there exists a subsequence such that u nk — )• u in C([0,T], L P (R )) for all 
1 < p < co. If additionally Uq G C°(M d ) and Uq — > uq locally uniformly, then we also have u nk — > u 
locally uniformly on [0,T] x R d . Moreover, T + {uq) = Tq < liminf n _> 00 T + (uq). 

Proof. From the local existence theory and the fact sup n ( 1 1 itg 1 1 cao + IKo Ik) < 00 > ^ ^ s assure d that 
Tq > (see [2] or [9]). Let < T < Tq. By the precompactness arguments of the local existence 
theory (see [21 H]) we may extract a subsequence {u nk } which converges in C([0, T]; L p (R d )) for all 
1 < p < oo to a weak solution, and by uniqueness of weak solutions, the limit must be u(t). In 
particular, we may extend u(t) to include any time interval such that 

liminf sup ||un(t)||oo < °°- 
te[o,T] 

By the proof of Lemma [5] and Theorem 6.1 of [21j . if Uq — > uq locally uniformly, it moreover 
follows that u nk (t) — >• u(t) locally uniformly up to extraction of an additional subsequence. By the 
continuation criterion in Theorem [TJ it is necessary that To < limin^^oo T + {uq) but there is no a 
priori reason for them to be comparable (for instance, To could be finite, but T + {uq) = oo). From 
the first part of the lemma, u(t) exists on all compact time intervals of [0, To). Since u nk {t) — > u(t) 
locally uniformly on this time interval, if To < oo, then we necessarily have liminf^ *Tq ||^(i)||oo = oo 
and therefore T + (uo) = To. This proves the lemma. □ 

3.2 Mass Comparison 

In this section we develop the comparison arguments and prove Theorem [8} 
First, note that u(t) is a solution to the transport equation 

d t u + V • (u-x) = 0. (3.7) 
H 

By construction, for all t, u(t,x) is also a weak solution to the scale-invariant problem (Remark [I]), 

— ^— V- (uVM*u) = Au m . (3.8) 
a(R ) 
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We have the following lemma which describes the PDE satisfied by the mass function corre- 
sponding to pit). 

Lemma 8 (Evolution of Mass Function). Let p(t,x) be a smooth radially symmetric solution to 
(3.3). Then, M(t,r) := J,, <r p(t,x)dx satisfies 

a t M it ,r> = (*^>Y" + ^%M W (3.9) 

\ ar a 1 / ayr) ar a 1 

where a is the surface area of the unit sphere in M. d . 
Proof. By radial symmetry 

p(t,r) = —^d r M(t,r) 
o~r a 1 

and by the divergence theorem and the radial symmetry of a, 

air) ! d r c(t,x)dS = — ^-Mit^r). (3.10) 
J\x\=r crr d 1 

Again using the divergence theorem and radial symmetry, 

8 t M(t, r)= f ~ • ( Vp m - p}- 2 / d pVc) dS 

J\x\=r \X\ V J 

= ar d-l dr ( 9rM (\ r) ) m - ^-d r M(t, r) [ d r c(t, r)dS 

V or d ~ v ) ar d 1 J\ x \ =r 

V ar a ~ L J ayr) ar a 1 

□ 



Suppose uq satisfies (1.10) and let p e (t) solve (3.3) with initial data 



e 

p € {0) = p + e 



-|x| 2 /4 

(Air) d / 2 ' 

By Lemma[6j p e (t,r) remain smooth and positive on their time interval of existence (0, T+(p e (0))). 

We are now ready to state and prove the mass comparison result which will complete the proof 
of Theorem [H 

Proposition 2. Suppose that T < min(T*, T+(p e (0))) and let M e (t,r) := J,, <r p e (t, x)dx . Further 
suppose M(0,r) < M t (0,r) for all r > 0. Then we have 

M(t, r) - M e (t, r) < in [0, T] x R d . 



Proof. Similar to the proof of Lemma pi by (3.8) we have, 



^(^r + ^^M W = 0, (3.11) 
V or 6 -' 1 J a(R ) ar d ~ x v ; ' y ' 
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and by (3.7), 



d t M(r) 



\x\=r 



u{r)-RdS 
R 



R 



Rd r M(r). 



(3.12) 



For notational simplicity, define M(t,r) := M e (t,r). 

Consider the space-time region (t,r) S Qt, where 

Q T = {(t,r):t£ [0,T], r 6 [0, R(t)}} . 

As M{u) = M(t,r) in \x\ > R(t), we need only prove the comparison result in Qt, from which the 
result on (0, T) x R d follows since M(t, r) is increasing in r. 

For a given constant A > (to be chosen later), let us consider the function 

f(t, r) := (M(t, r) - M(t, r))e~ A * in [0, T] x R+. 

Note that /(0,r) < 0. If f(t,r) < in Qt there is nothing to prove, so suppose that it is positive 
somewhere. Then f(t,r) has a strictly positive maximum in Qt, which is achieved at some point 
(t+, r*). Necessarily, > 0. If = R(t+) then by r* being the location of the maximum, we must 
have 

= d r M(i*,r*) > r M(t*,r*) = <rrf-V(t*,r*). 
Since ii e is strictly positive, it follows that r* < R(t+). This implies that due to maximization, 

(A) d t (M(t^ r*) - M(t*, r*)) > A(M(i*, r*) - M(i*, r*)). 

(B) 9 r M(t*,r*) = ^.M(t*,r*). 

(C) d rr M(U,n) > d rr M(U,r*). 



Using the mass equations satisfied by each function we have, by (3.9), (3.12) and (3.11), 

~~ "' M 



d t {M-M){U,r, 



urt X d, 



d r M 



or+ 



.d-l 



+ 



1 



.d-l 



or. 



.d-l 



Or 



d r M 



or 



d-l 



cl(Ro) or. 
p~ 2 l d M 



a ( r *) 



or 



d-l 



d r M 



d r M - d r M 



R(U) 



R(Q. 



where all terms are evaluated at (i*, r+). By (B) and (C), we can order the higher order nonlinearity 
coming from the diffusion as well as relate the advection terms, 



d t (M - M)(t*,r*) < 



d r M 



a(R 



1 M " 



-2/d 



M 



d r M 



-R(U 



Using that r* < Rq and by assumption a(r) is non-decreasing on r 6 [0, Rq] we have, 



3 t {M - M){U,r*) < 



1 d r M 



+ d r M 



(1 - /i" 2 / d ) M r* 
W*) art 1 ~ R(U) 



(i*,r*). 
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Since V is radial and non-increasing, we have 

d r (r 1 - d d r M(t,r)) = a- l d r (u{t,r)) < 0. 

In addition we have M(t,0) = and M(t,R(t)) = M c . Note that h(r) := M c (^) d solves 
drir^drhir)) = in (0,R(t)) with boundary data h(0) = and h(R(t)) = M c . Therefore it 
follows from the weak elliptic maximum principle, 



M(t,r) > M c 



R(t) 



Using the above observation we have 

d t {M - M){U,r+) < -i ^ ^- (M-M)(4,r, 



+ 



d-l 



[1 - pT 2 / d ) M M 



Due to (3.5), and using the fact that a(r) is non-decreasing, we have 
d t (M-M)(U,r+) < 1 dr ™ (M-M)(U, rir ) 
+ 8 r M ^ - .J (M-M) 
Since A(M — M) < d t (M — M), the result follows by choosing 



A > sup \\u 

te[o,T] 



1 



t(0)/i 2 /*' 



We may now prove Theorem [8) 



(3.13) 



□ 



Proof. (Theorem [8]) Let uq satisfy the hypotheses of Proposition [T] and let {p € } e>0 and T* be given 
as above. By the hypotheses of Proposition [TJ 

M(0,r) < M(0,r) < M e (0,r) for r G [0,oo) for any e > 0. 

Let To be defined as in Lemma[7J which satisfies T + (po) = To < liminf e _s.o T + (p e (0)). Therefore, 
it suffices to show To < T*. 

To this end, suppose To > T*, which implies that {p e (t)} and p(t) exist on [0,TJ for sufficiently 
small e. Moreover, by Lemma [7J there exists a sequence p tk — > p in C([0, TJ; L 1 (M d )) and locally 
uniformly. Combined with Proposition [2j this implies 



M(t,r) < M(t,r), 



(3.14) 



for all < t < XJ Since u concentrates at time TJ (3.14) implies that p(t) must also concentrate 
at T*, contradicting the assumption T + (uq) = Tq > T*. □ 
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We briefly sketch a proof of Theorem [9| 



Proof. (Theorem [9]) Let uq be as in the statement of Theorem [9j We only need to verify (1.10) in 
order to apply Proposition jlj Let R\ be such that Mq = f\ x \<R 1 'Uodx > M c . Then (1.10) holds 



for any r > R\. By assumption, M C M 1 uq is strictly positive on some compact ball {\x\ < r%}. 
Hence we may choose Rq sufficiently large such that Rq c1 V(Rq 1 x) < M c Mq 1 Uq(x) for \x\ < r\ and 



therefore (1.10) holds up to at least r = r\. As for r\ < r < R\, M c Mq f,, <r uo(x)dx is non- 
decreasing and hence bounded below on the compact annulus by the value at r = n. Hence, we may 
choose Rq even larger to ensure that ( |1.10 ) holds also for n < r < R\ and therefore everywhere. 



Hence, we may apply Proposition [T] and the result follows. 

We now prove that we may construct blow up solutions with arbitrarily large initial free energy. 
We follow a similar procedure as Lemma 3.7 in [9]. Let Vm C L\ n L rn be the set of non- negative, 
radially symmetric non-increasing functions in L 1 n L m with mass M. By the above reasoning, if 
uq G Vm is continuous with finite second moment then the associated solution u(t) to the scale- 
invariant problem with u(0) = uq blows up in finite time. Now we prove that 



Suppose for contradiction that 



sup F{h) = +00. 

h&V M 



A := sup F{h) < 00. (3.15) 

/eV M 



Following the same scaling argument as in Lemma 3.7 of [9j, we may use the HLS (1.6) to show 



(3.15) implies a reverse Holder- type inequality for any h G L\ n L m which is radially symmetric 



non-mcreasmg 



2 l d < \\h\\ 2 
>.M \\tl\\ 



lmll'"lll r^M \\ n \\2d/(d+2)- 

However, this inequality is definitely false, as m > 2d/(d + 2) implies we may easily construct a 
sequence of functions with uniformly bounded L 2d /( d + 2 ) norm and unbounded L m norm. Indeed, 
consider f s = (5+\x\)- a l B{0tl) (x) with a G (d/m, (d+2)/2) with 5 G [0, 1]. Then \\f s \\i > ||/i||i > 
but ||/s||2d/(d+2) is uniformly bounded and lim^^o ll/<s||m = 00. Hence it follows by contradiction 
that A = +00. By density, we may restrict to continuous functions with finite second moment in Vm 
and show that there are solutions to the scale-invariant problem with initial free energy arbitrarily 
large which blow up in finite time. □ 



4 Global existence at critical mass 

In the scale-invariant case, the free energy no longer controls the entropy if M(u) = M c . Indeed, 
consider a family of extremals (from Theorem p5J), Vr = R~ d V{R~ l x) with R — >■ 0, which concen- 
trates into a Dirac mass but all satisfy -F(Vr) = 0. Hence unlike the subcritical case, a proof of 
global existence at the critical mass requires more information than what is provided by the energy 
dissipation inequality alone. 

In the homogeneous case, the mass comparison principles with a supersolution prove global 
existence of radial threshold solutions |45] , a result which extends to general threshold solutions by 



symmetrization inequalities |19 | 120 } [30]. However, it is important to note that (1.1) does not satisfy 
the same comparison principles as the homogeneous problem and in particular, radial solutions 
which are initially monotone are not guaranteed to remain so. Among other things, this implies 
that a mass supersolution alone cannot be used to imply global existence. However, Theorem 
[6] provides a kind of rigidity to threshold blow-up solutions, as all of the mass in the solution 
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must concentrate into a single point where a(x) achieves the minimal value. It is this additional 
information which allows us to prove Theorem [7j as now we only need to rule out a very specific 
behavior. Indeed, we use a mass supersolution to show that any radially symmetric threshold 
solution cannot concentrate into the origin, even if the blow-up is in infinite time. Therefore, by 
Theorem [6j all radially symmetric critical mass solutions must remain uniformly bounded. 

Certain details of the mass comparison argument used to prove Theorem [7J differ from those 
used to prove Theorem [8j In fact, the proof of Theorem [7j is significantly easier. In this section 
we use the extremals as mass supersolutions, as opposed to the previous section where our barriers 
were subsolutions. Define 



min a(x), 
x&R d 



and our supersolution to be the extremal 



HT V (S)" (41) 

where now R > is a sufficiently small constant which will be fixed later in the proof. Note as well 
that M{u) = M c . As above in we define the mass distributions 



M(t,r)= / u{t,x)dx, M(t,r)= I u(t,x)dx. 

J\x\<r J\x\<r 

In contrast to R is chosen such that M(0, r) < M(r, 0) for all r € [0, R], which is always possible 
since no is continuous. We now show the following proposition. Since the proof is similar to the 
proof of Proposition [2] we only sketch it. In light of Theorem [6j this proposition implies Theorem 
[JJ as it rules out the finite or infinite time concentration of mass. 

Proposition 3. Suppose that T < min(T + (uo)). Further suppose M(0, r) < M(Q,r) for all r > 0. 
Then we have 

M(t,r) - M(t,r) < in [0,T] x R d . (4.2) 
Proof. As in the previous section, we must regularize the initial data. Hence, fix e > 0, define 

1 f \x~y\ 2 

u e (x) = — e 4, u (y)dy. 

Note that we use a different regularization here than in the proof of Proposition[2]in order to preserve 
the mass. Denote the associated weak solutions u e (t), which are smooth and strictly positive until 
blow-up time T+(mq) by Lemma [6j We will deduce (4.2) independent of e and hence Lemma [7] shows 
that the inequality holds also for u(t). 



The proof that (4.2 ) is satisfied independent of e is a simple variation on the proof of Proposition 
[2j made significantly more straightforward by the simpler barrier u(x, t) which here does not depend 
on time. The difference is that u(x, t) plays the role of supersolution here, and so we are looking for 
the positive maximum of (M(t,r)-M(t,r))e- Xt . The main place the argument differs significantly 
from Proposition [2] is along the contact line \x\ = R at the edge of the support of u. We can 
rule out that a positive maximum occurs here for the simple reason that the barrier and the weak 
solutions have the same mass. Hence the solutions u e cannot possibly have strictly more mass in a 
ball of I a; I < R than the barrier. Therefore, any positive maximum must occur somewhere in the set 
< \x\ < R. From this point, we may continue as in the proof of Proposition [2j Note that unlike 
the proof of Proposition [2j the presence of "f(x) does not pose a problem here. □ 
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Remark 10. Theorem [6] shows that any solution with critical mass can only blow up by con- 
centrating all of the mass into a point where a(x) achieves the minimum value. Hence with the 
additional constraint of radial symmetry, Theorem [6] directly implies Theorem [7] without the need 
for Proposition [3] unless mina(x) = a(0). 
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Appendix: Pressure Form Comparison 

By Lemma [5j for any e > we have 

L := sup |Ac| < 00. (4.3) 

e<t<T*-e 

Let us define the pressure form of u: 

v = -^—(uT-\ (4.4) 
m — 1 

Then formally v solves the following equation: 

vt = (m — l)vAv + \Dv \ 2 + X7v ■ Vc + (m — l)vAc. 

We proceed to prove the "viscosity solution" property of the pressure v. The notion of viscosity 
solutions are first introduced for Hamilton- Jacobi equations by Crandall-Lions ([H]) and later for 
fully nonlinear elliptic-parabolic equations (PZ]) as wen as f ree boundary type problems (see [T3] . 
|14j and |29] for porous medium-type problems.) The advantage of the approach lies in pointwise 
control of solutions and, in our setting, their free boundaries. More specifically we will show that 
the initially positive solutions cannot touch down to zero at later times, i.e. that contact lines 
cannot be nucleated. 

Since c is not C 2 up to the zero set of u, Ac is not well-defined on the free boundary 
d{u > 0} = d{v > 0}. This causes a technical problem for directly applying a standard notion of 
viscosity solutions to v. Hence we will directly prove the necessary properties to be used in our 
analysis in the next section. 

Definition 3. For nonnegative functions u and v defined in a small neighborhood £ of (xo,to), we 
say 

(a) u crosses v from below at (xo,to) if 

u < v in S n {t < to} and u(xq, to) = v(xq, to). 
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(b) u crosses v from above at (xq, to) if 

u > v in E n {t < to} and u(xo, to) = v(xq, to). 



Proposition 4. For any given domain E C M. d x [0,T* — e), Ze£ be a nonnegative continuous 
function in E which is C 2,1 in {cp > 0} ; u>i£/i \D<j)\ > on 9{i^> > 0}. Let v be the pressure form of 
u as defined in (4.4). Then the following holds: 



(a) Suppose v crosses (f> from below at (xo,to) in {u > 0} n {t < to} in E. Then we have 

4>t-{m- l)4>A(j) - \D(/)\ 2 - V<f> ■ Vc - (m - l)Lcj> < 



(b) Suppose v crosses 4> from above at (xo,to) in {u > 0} n {t < to} in E. Then we have 

(j> t -(m- 1)0A0 - \D(f)\ 2 - V0 • Vc + (m - > 0. 



i^ere i/ie constant L is as given in (4.3) 



Proof. 1. Note that it and t> are smooth in their positive set, and there the result follows easily. 
Hence, the only difficult case is when (xo,io) 6 d{v > 0}. 

2. Let us take c e := c* r/ t , where r/ t is the standard mollifier. Let u e be the weak solution of 

(u t ) t = A(u e ) m + V-(u e Vc e ), 

with initial data u e (x, 0) = u(x, 0). Since c e is C 2 , it follows from [29] that n e is a viscosity solution in 
the sense defined in therein. In particular, it is shown in [29] that the statements in the proposition 
hold for v t : the pressure form of u e . Below we will approximate u by u e to prove the proposition. 
Since c e is uniformly bounded in C 1 ' 1 norm, u e is equi-continuous due to Theorem 6.1 of |21j . Using 
this fact, parallel arguments leading to Proposition 3.3 of [30J yields that u e uniformly converges to 
u, and thus v e to v. 

3. Let us now show (a) when v crosses a nonnegative function <ft 6 C 2,1 ({(j) > 0}) from below 
at (xo,to) G d{v > 0}. Let us perturb </> so that v really crosses 4>, not just touching. This can be 
done by replacing eft with 

4>(x, t) = (4>(x, t) - a(t -t - &))+, 

where a and b are small positive constants. 

4. If (a) fails then, since ^(cco^o) = 0, then (j> satisfies 

Ot - (m - l)(pA(p - \D(p\ 2 - V0 • Vc - (m - l)L0)(x o , to) > 0. 

Now let us pick a small 5 > and take 4>s{-,t) := (0 + )(-, i) * 7/5 + m(o~) where n(x) is a standard 
mollifier which is smooth and has exponential decay at infinity, and m(5) is a constant. Choose 
m(5) accordingly so that v is strictly below <f>s at t = to but crosses 4>$ from below at t = to + 0(5). 
Note that this is possible because rj does not have a compact support. 
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Then due to continuity of the derivatives of <j> in its support and the corresponding convergence 
of 4>s to (j), <ps satisfies 

((f>s) t - (m - l)<t> S A<t>5 ~ \D(f>s\ 2 - Vfa ■ Vc e - (m - l)L<j) & > 0. (4.5) 

in O((5o)-neighborhood of (xo>*o) if e 3 <5 << <5o- 

Since u e converges uniformly to ti as e -> 0, u e crosses cfo from below at (xs,t$), which lies in 
O((5o)-neighborhood of (xo>*o) if e and <5 are chosen small enough. 
Note that at (xq, to) we have 

(v e )t > {D(f> s )t, \Dv £ \ = \D(f>s\ and Av e < Afa. 



this contradicts (4.5) and the fact that v £ satisfies 

(v e ) t - (m - l)v f Av € - \Dv e \ 2 - Vv e ■ Vc e - (m - l)Lv e < 

□ 

Corollary 1 (Local comparison in pressure variable). In any given parabolic, cylindrical neighbor- 
hood E, Zef 4> be a C 2,1 function in {4> > 0} io#/i > on d{4> > 0}. 

(a) Suppose that <fi satisfies 

4> t -{m- l)cf>A(t> - \D<p\ 2 - Vcf> • Vc - (m - l)L<j) > in E. 

T/ien u cannot cross <p from below in E. 

f&j Suppose (j) satisfies 

4> t -(m- l)<f)A<j> - \D4>\ 2 - V4> • Vc + (m - 1)L0 < in E. 

Then v cannot cross 4> from above in E. 

An immediate consequence of the above proposition is the preservation of positivity for u. 

Lemma 9. Let u(x, t) be the weak solution associated with the strictly positive continuous initial 
data uq(x) > 0. Then u(x,t) is strictly positive everywhere up to the blow-up time T + (uq). 

Proof. Let v be the corresponding pressure form of u. Let us recall that the Barenblatt profile is 
given as 

B{x,t) :=r* (c-k^ 
where C > is a positive constant and 

d(m - 1) _ A _ A 
~ d(m- 1) + 2'^ ~ d' ~ 2d' 

B{x, t) then solves the porous medium equation in its pressure form in the viscosity sense (see e.g. 
|29]): 

B t -(m- 1)BAB - \DB\ 2 = 0. 

Let us now define 

B(x, t) := e~ Mt sup B(x, t) for < t < 1. 

yeB M _ Mt {x) 
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Then due to Proposition 2.13 in [29J B satisfies 

B t -(m- 1)2? AS - \DB\ 2 + M\DB\ + MB < 
for < t < 1. Let us choose 

M = (m- l)max(||Vc|| L oo,L) . 
(Note that the first term in above upper bound is bounded before the blow-up time). 

Since B(x,t) vanishes uniformly to zero as C — > 0, so does B. Hence for any r > one can 
choose C = C(r) sufficiently small so that B(x,t) < uo- Then Corollary [l] yields that 

B(x, t + t) < u(x, t) for < t < 1. 

Since r can be arbitrarily large and B has its support expanding to the whole domain as r grows 
to infinity, we conclude that u is strictly positive for < t < 1. 

We can iterate above argument up to the blow-up time to conclude (since the solution and L 
remain bounded until blow-up). 

□ 
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